MICROLOCAL ASPECTS OF BISTATIC SYNTHETIC 
APERTURE RADAR IMAGING 



Abstract. In this article, we analyze the microlocal properties of the lin- 
earized forward scattering operator F and the reconstruction operator F*F 
appearing in bistatic synthetic aperture radar imaging. In our model, the 
radar source and detector travel along a line a fixed distance apart. We show 
that F is a Fourier integral operator, and we give the mapping properties of the 
projections from the canonical relation of F, showing that the right projection 
is a blow-down and the left projection is a fold. We then show that F*F is a 
singular FIO belonging to the class I 3 ' . 
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1. Introduction 

In synthetic aperture radar (SAR) imaging, a region of interest on the surface of 
the earth is illuminated by electromagnetic waves from a moving airborne platform. 
The goal is to reconstruct an image of the region based on the measurement of 
scattered waves. For an in-depth treatment of SAR imaging, we refer the reader 
to [3l [2] ■ SAR imaging is similar to other imaging problems such as Sonar where 
acoustic waves are used to reconstruct the shape of objects on the ocean floor 

EG3Hg|. 

In monostatic SAR, the source and the receiver are located on the same moving 
airborne platform. In bistatic SAR, the source and the receiver are on independently 
moving airborne platforms. There are several advantages to considering such data 
acquisition geometries. The receivers, compared to the transmitters, are passive 
and hence are more difficult to detect. Hence, by separating their locations, the 
receivers alone can be in an unsafe environment, while the transmitters are in a safe 
environment. Furthermore, bistatic SAR systems are more resistant to electronic 
countcrmcasures such as target shaping to reduce scattering in the direction of 
incident waves [17] . 

In this paper, we consider a bistatic SAR system where the antennas have poor 
directivity and hence the beams do not focus on targets on the ground. We assume 
that the transmitter and receiver traverse a 1-dimensional curve and the back- 
scattered data is measured at each point on this curve for a certain period of time. 
As in the monostatic SAR case [21] , with a weak scattering assumption, the linear 
scattering operator that relates the unknown function that models the object on the 
ground to the data at the receiver (see [35]) is a Fourier integral operator (FIO) [TBJ 
[IHlIll [24]. Now if T is an FIO, the canonical relation Ajr associated to J- tells us how 
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the singularities of the object are propagated to the data. The canonical relation 
Ajf* of the L 2 adjoint J-* of J- gives us information as to how the singularities in 
the data are propagated back to the reconstructed object. The microlocal analysis 
of singularities of the object is then done by analyzing the composition A.f» o Ajr. 

Such an analysis for monostatic SAR has been done by several authors [HI [SJ 
[7] and is fairly well understood. In their work [22], Nolan and Cheney showed 
that the composition of the linearized scattering operator with its L 2 adjoint is 
a singular pseudodiffcrential operator (^DO) and it belongs to class of Fourier 
integral operators with two cleanly intersecting Lagrangians. Felea in her works 
[5[ [7] further analyzed the properties of the composition of these operators. 

In this paper, we do a similar analysis for the bistatic SAR imaging problem. 
Given the complications that arise in treating arbitrary transmitter and receiver 
trajectories, in this paper, we focus on the case where the transmitter and receiver 
are at the same height above the ground, traverse the same linear trajectory at 
the same constant speed and spaced apart from each other by a constant distance. 
Furthermore, we assume that the object to be imaged is on the ground, which for 
simplicity, we will assume is flat. Since the measured data is two-dimensional, it is 
reasonable to expect that we can reconstruct a two-dimensional object. 

The outline of the paper is as follows. Section [2] focuses on the preliminaries. 
Here we give the linearized scattering model for bistatic SAR and definitions of 
singularities and distributions that belong to I p ' 1 classes and important results on 
distributions belonging to this class that are required in this paper. In Section [3j 
we undertake a detailed study of the canonical relation associated this FIO. This 
is the content of Theorem 13.51 Then we study the reconstruction operator, that 
is, the composition of the bistatic scattering operator with its adjoint, and show in 
Theorem 14.31 that this operator belongs to the class I 3,0 . Our proof of 14.31 follows 
the ideas of [5J Theorem 1.6]. Several identities required to prove Theorem 14.31 arc 
provided in the Appendix. 



2.1. The bistatic linearized scattering model. We assume that a bistatic SAR 
system is involved in imaging a scene. Let Jt(s) and jr(s) for s S (soj s i) be the 
trajectories of the transmitter and receiver respectively The transmitter transmits 
electromagnetic waves that scatter off the target, which are then measured at the 
receiver. We are interested in obtaining a linearized model for this scattered signal. 

The propagation of electromagnetic waves can be described by the scalar wave 
equation: 



where c is the speed of electromagnetic waves in the medium, E(x, t) is the electric 
field and P(t) is the transmit waveform sent to the transmitter antenna located at 
position 7t(s)- The wave speed c is spatially varying due to inhomogencities present 
in the medium. We assume that the background in which the electromagnetic waves 
propagate is free space. Therefore c can be expressed as: 
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where the constant cq is the speed of light in free space and V{x) is the perturba- 
tion due to deviation from the background, which we would like to recover from 
backscattered waves. 

Since the incident electromagnetic waves in typical radar frequencies attenuate 
rapidly as they penetrate the ground, we assume that V(x) varies only on a 2- 
dimcnsional surface. Therefore, we represent V as a function of the form 

V(x) = V{x)5 (x 3 ) 

where we assume for simplicity that the earth's surface is flat, represented by the 
x = (x±,x2) plane. 

The background Green's function g is then given by the solution to the following 
equation: 

A - \d?) g(x, 7 t(s), t) = 5(x - M*))8(t)- 
We can explicitly write g as 

ff(*,^(.),*) = ^P^). 

4tt|x - 7t(s)| 

Now the incident field E ln due to the source s(x, t) = P(t)5(x — 7t(s)) is 



E m (x,t) = j g(x,y,t~T)s(y,T)dydT 

_ P(t-\x- lT ( S )\/c ) 
4ir\x — 7r(s)| 

Let E denote the total field of the medium, E = E ln + E sc . Then the scattered 
field can be written using the Lippman-Schwinger equation: 

(1) E sc {z,t) = J g{z,x,t-T)d?E(x,T)T(x)dxdT. 

We linearize ([T]) by the first Born approximation and write the linearized scattered 
wave-field at receiver location 7r(s): 

£fi c „(7fl(s), t) = j g( lR (s), x, t - r)d?E in (x, r)V (x)dxdr 

r 2 ) = f ^( t - T -\ x -lR(s)\/c ) f_ iu{T _ lx _ jT{a)l/co) uj 2 P (uj) 

J 47t|x-7_r(s)| V 47r|a:-7T(s)| 

x V(x)dLudxdr, 

where p is the Fourier transform of P. 

Now, integrating ([2]) with respect to r, a linearized model for the scattered signal 
is as follows: 

(3) d(s,t) := Ef^ R (s),t) = ( e-^-^*)) A {s,x^)V{x)dxdu, 



where 

(4) R(s,x) = \~f T {s)-x\ + \x-~( R (s)\ 
and 

(5) A(a,x,u) =p{uj)({Ati) 2 \ 1t ( s ) - x\\ lR {s) - x\)-\ 

This function includes terms that take into account the transmitter and receiver 
antenna beam patterns, the transmitted waveform and geometric spreading factors. 
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We will show in Section [3] in one important case, that the transform T that 
maps V to ((3|) is a Fourier integral operator associated to a canonical relation A 
(Proposition 13. lj) . and we will prove the mapping properties of A (Propositions 13.21 
and 13. 4j) . These mapping properties tell what T does to singularities. We now 
define the mapping properties we need. 

2.2. Singularities and I p > 1 classes. Here we give the definitions of the singular- 
ities associated with our operator T and its canonical relation (fTTj). and a class of 
distributions required for the analysis of the composition T with its L 2 adjoint. 

Definition 2.1. [TT] Let M and TV be manifolds of dimension n and let / : N — > AI 

be C°°. 

(1) / is a Whitney fold if near each m S M, f is either a local diffcomorphism 
or / drops rank simply by one at mo so that L = {m G M : rankdf = n— 1} 
is a smooth hypcrsurface at mo and kcrc?/(mo) <£. T mo L. 

(2) / is a blow-down along a smooth hypcrsurface L C M if / is a local dif- 
fcomorphism away from L and / drops rank simply by one at L, where 
the Hessian of / is equal to zero and kerdf C T(L), so that f\ L has one- 
dimensional fibers. 

We now define I p ' 1 classes. They were first introduced by Melrose and Uhlmann, 
[20] Guillemin and Uhlmann [15] and Greenleaf and Uhlmann [10] and they were 
used in the context of radar imaging in [22j [6j [7] . 

Definition 2.2. Two submanifolds M and TV intersect cleanly if Mfl N is a smooth 
submanifold and T(M D N) = TM n TN. 

Consider two spaces X and Y and let Ao and Ai and Ao and Ai be Lagrangian 
submanifolds of the product space T*X x T*Y. If they intersect cleanly (A , A x ) 
and (Ao,Ai) arc equivalent in the sense that there is, microlocally, a canonical 
transformation x which maps (Ao, Ai) into (Ao, Ai). This leads us to the following 
model case. 

Example 2.3. Let A = A t *r™ = {(x, ^; x, g)\x Si?", £, e R n \ 0} be the diagonal 
mT*R n x T*R n and let A x = {(x 1 , x n , 0; x', y n , 0)\x' e i?™^ 1 , f £ J^-^O}. 
Then, Aq intersects Ai cleanly in codimension 1. 

Now we define the class of product-type symbols S p ' l (m,n, k). 

Definition 2.4. S p - l (m, n, k) is the set of all functions a(z, £, a) G C°°(R m x R™ x 
R fc ) such that for every K C R m and every a £ Z™ , f3 £ Z" , 7 € Z^ there is CK, a ,/3 
such that 

|^d>(z,£,a)| < CM (1 + |?ir l/3| (l + |a|)H-vl )V ( z ^, r ) e K x R" x R fe . 

Since any two sets of cleanly intersecting Lagrangians are equivalent, we first 
define I pJ classes for the case in Example 12.31 

Definition 2.5. [T3] Let I p,1 (Aq, Ai) be the set of all distributions u such that 
u = u\ + U2 with Mi e Cq 30 and 

with a e where p' =p — % + \ and V = l-\. 
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This allows us to define the J p ''(Ao,Ai) class for any two cleanly intersecting 
Lagrangians in codimension 1 using the microlocal equivalence with the case in 
Example [231 

Definition 2.6. [T5] Let I p > (Aq,Ax) be the set of all distributions u such that 
u — u 1 + u 2 + J2 v t where ui £ I p+l (A \ Ax), u 2 G / p (Ai \ A ), the sum J2 v i 
is locally finite and Vi = Fwi where F is a zero order FIO associated to Y -1 , the 
canonical transformation from above, and Wi G I p ' (Aq, Ai). 

This class of distributions is invariant under FIOs associated to canonical trans- 
formations which map the pair (Ao, Ai) to itself. By definition, F G I P ' Z (A , Ai) if its 
Schwartz kernel belongs to I p ' l (A , Ai). If F G I P ''(A , Ai) then F G I p+l (A \ A x ) 
and F G 7 p (Ai \ A ) [TS]. Here by F G I P+ '(A \ Ai), we mean that the Schwartz 
kernel of F belongs to 7 p+/ (Ao \ Ai) microlocally away from Ai. 

3. Transmitter and receiver in a linear trajectory 

Henceforth, let us assume that the trajectory of the transmitter is 

7t : (s ,si) ^M 3 , j T (s) = (s + a,0,h) 

and the trajectory of the receiver is 

7i?,(s) : (s , Si) -> M 3 , 7 fl (s) = (s - a, 0, /i). 

Here a > and ft, > arc fixed. From Equation ©, the linearized model for the 
data at the receiver, for s G (so,si) and i G (io,ii) is 

(6) d(a,t) = / e- iti, ( t -* (|x - 7T(s)l+|x - 7RWI) )^( S)a; ,a;)V r ( : r)d a; da;. 



We multiply d(s,t) by a smooth (infinitely differentiable) function f(s,t) supported 
in a compact subset of (so, si) x (^o, This compensates for the discontinuities in 
the measurements at the end points of the rectangle (sq, s\) x (to, t\). For simplicity, 
let us denote the function / • d as d again. We then have 

(7) d(s,t)= f e-^ {t -^ R{s ' x)) A(s,t,x,u)V(x)dxdLo, 



where now A(s,t, x, oj) = f(s,i)A(s,x,ui). 

Our method cannot image the point on the object that is "directly underneath" 
the transmitter and the receiver. This is, if the transmitter and receiver are at 
locations (s + a, 0, h) and (s — a,0, h), then we cannot image the point (s,0, 0). 
Therefore we modify d in Equation ([7]) by multiplying by another smooth function 
g(s, t) such that 



g = in a small neighborhood of < s, 2 : sq < s < Si >. 

H c ° / J 

For simplicity, again denote g ■ d as d and g • A as A. Consider, 
(8) FV(s, t) := d(«, t) = / e-^-^-y^l+^-^i^Ais, t, x, oj)V{x)dxdu. 



For simplicity, let us denote the (s,t) space as Y. 

We assume that the amplitude function A satisfies the following estimate: For 
every compact K G Y x X and for every non-negative integer a and for every 
2-indexes /? = (/?i,/?2) and 7, there is a constant C such that 

(9) \dZdl-dtdZA{s,t,x,u))\ < c(i + M) 2 - tt . 
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The phase function of the operator F, 
(10) ip(s,t,x,(j) = -u(t- — {\x - 7T (s)| + \x - j R (s)\) 



CO 

is homogeneous of degree 1 in w. 

We now analyze some properties of the canonical relation of the operator F. 

Proposition 3.1. F is a Fourier integral operator of order 3/2 associated with 
canonical relation 

(11) 



c o VI^-TtIs)! \x-j r (s) 



f uj I x\ — s — a x\ — s — a \ uj X2 X2 

\xi,x 2 , 



co\\x-f T (s)\ \x-j R (s)\J Co \\x-j T (s)\ \x-j r (s)\ 
■■Cot = \j {x\ - s - a) 2 + x\ + h 2 + (xi -s + a) 2 +x\ + h 2 , w^oj. 
Furthermore {x\,X2, s,uj) is a global parameterization for A. 

Proof. This is a straightforward application of the theory of FIO. Since ip in (| 10[) is 
a nondegcnerate phase function with d x ip and d s .ti}) nowhere zero and the amplitude 
A in ([H) is of order 2, F is an FIO [16]. Since the amplitude is of order 2, the order 
of the FIO is 3/2 by [HI Definition 3.2.2]. By definition [HI Equation (3.1.2)] 

A = {(s,t, d Stt ip(x,s, t,u)), (x, -d x ip(x, s,t)) : d u ip(x,s,t,u) = 0}. 

A calculation using this definition establishes ((TTj) . Finally, it is easy to see that 
{x\,X2, s,uj) is a global parameterization of A. 

□ 

In order to understand the microlocal mapping properties of F and F*F, we 
consider the projections tt l : T*Y x T*X T*Y and ir R : T*Y x T*X T*X. 

Proposition 3.2. The projection ttl restricted to A has a fold singularity on £ :~ 
{(x,0,s,u) : uj ^ 0}. 

Proof. The projection 7r^ is given by 

(12) TTl(xi,X2,S,L)) = 

I I/, r \ i i / \ i \ u) ( X\ — s — a xi — s + a 

= Is, — (\x - 7 T (s)| + \x- Jr{s)\) , 



co c VN-7t(s)| \x-j r (s) 

We have 

/ 1 \ 



(ttx). = 



Then 



1 / xi — s — a. I x\ — s-\-a 

c I \x--y T (s)\ \x-j R (s)\ J c \\x-'y T (s)\ 1 |:e-7.r(-s) 



c ^ |s- 







u / X2+I1 2 x?,+h 2 \ / (zi— s — 0)3:2 i (ii-s+a)i2 

"co VI^-7T(s)| :i + | K - 7H ( S )|^ co V k-7r(s)| 3 |x- 7R (s)| 3 

y °-iy 



det(7r ) - —x ( 1 + 1 \ A + (xi ~s) 2 +^ + fe 2 -a 2 
L * c § 2 V|x - 7t(s)| 2 |x-7 fl (s)| 2 / V |a;-7T(s)||a;-7fl(s)l 
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The proposition then follows as a consequence of the following lemma. 
Lemma 3.3. The term 

1 (xi - s) 2 + x\ + h 2 - a 2 
\x - j T {s)\\x - j R (s)\ 

is positive for all i£l 2 , s G R and h and a positive. 

Proof. This is a straightforward calculation that is made simpler if one lets 5* = 
X\ — s, T = \/ x\ + h 2 then puts the term over a common denominator. Finally, one 
shows that the numerator is positive by isolating the square roots (absolute values) , 
squaring, and simplifying to infer that, since AT 2 a 2 > 0, the lemma is true. □ 



Now returning to the proof of the Proposition 13.21 we have that det(7Ti) H < = if 
and only if xi = 0. Hence det(7Ti)* vanishes on the set E and Lemma 13.31 again 
shows that d(dct(7T£)*) on E is non- vanishing. This implies that ttl drops rank 
simply by one on E. 

Now it remains to show that TE n Kernel^^)* = {0}. But this follows from the 
fact that Kernel(TTi)* = span(g|-), but TE = span(g^-, J^, -^). This concludes 
the proof of the proposition. □ 



Proposition 3.4. Consider the projection ttr : T*Y x T*X — s- T*X. The restric- 
tion of the projection to A has a blowdown singularity on E. 

Proof. We have 

(13) ttr(xi,X2,s,lj) = 



uj / x\ — s — a X\ — s — a 

•^1 *^2 I 

' ' c V>-7t(s)| \x-~/r(s)\ 



Co 



•'2 



•''2 



\x-~f T {s)\ \x-j R (s) 



Now 



/I 








1 






V 


* 



/ x%+h 2 
\ \x-7 T (s] 
(xi—s — a)x 



)l 3 
a)x 2 _ 
c \ \x-~/ T (s)\ 3 



x~+h' 



{xi~s+a)x 2 




\x—y R (s)\ J J 



From this we see that Kernel(7Tij)* C TE. Since det((7Tfl)*) = det((7T£)„), ttr drops 
rank simply by one along E. Therefore the projection ttr has a blowdown singularity 
along E. □ 



We summarize what we have proved in this section by the following theorem: 

Theorem 3.5. The operator F defined in ([8]) is a Fourier integral operator of 
order 3/2. The canonical relation A associated to F defined in (jlip satisfies the 
following: The projections ttl and ttr defined in (fT2"|) and (fT3"|) are a fold and 
blowdown respectively. 
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4. Image Reconstruction 

Next, we study the composition of F with F*. This composition is given as 
follows: 



x A(x, s, t, uj)A{y 1 s, t, Lo)dsdtdujdujdy 

After an application of the method of stationary phase, we can write the kernel of 
the operator F*F as 

K( X , y)= f e i #(l«'-T*(»)l + l»-'YHW|-(|x-7 ; rWI + |x-7i l WI))^( a . ) y; Sj w ) dsdw . 

Therefore the phase function of the kernel K(x, y) is 
uj 

(14) </>(x,y,s,u) = — (\y-rr(s)\ + \y- 1r{s)\ - (f~7t(s) + \x- 7a(s)|)) . 
co 

Proposition 4.1. 

WF(K)' c A U A 

w/iere A := {(xi, x 2 , £i, 6; x 2 , £i, 6)} arid A := {(xi, x 2 , £i, £2; xi, —%2, £1, -&)}- 
_ffere /or a point x = (xi,X2), £/ie covectors (£i,£ 2 ) ore non-zero multiples of the 
vector (—d Xl R(s,x),—d X2 R(s,x)), where R is defined in Q. 

Proof. Using the Hormandcr-Sato Lemma, we have 

WF{K)' C 

ui I xi — s — a x\ — s + a \ w / x 2 x 2 

C < I Xi,X 2 , 



c \\x-jt{s)\ \x-jr(s)\J' co \\x - tt(s)| |x-7r(s)| 

w/?/i-s-a j/i-s + a\ cj / y 2 2/2 
2/1,2/2, 1 



c o VI2/-7t(s)| |y-7fl(s)l/ c vIj/ _ 7t(s)| \y-~/R(s)\ 
|x - 7r(s)| + |x - 7h(s)| = |y - y T (s)\ + \y- 1r(s)\, 
xi — s — a x\ — s + a yi — s — a yi — s + < 



|x-7t(s)| |x-7r(s)| |y-7r(s)| l2/-7/?(s)| : 



w 7^ 



We now obtain a relation between (xi, x 2 ) and (j/i, y 2 ). This is given by the following 
lemma. 

Lemma 4.2. For al/ s, i/ie set of all (xi,X2), (2/1,2/2) satisfy 

(15) |x - 7t(s)| + |x - 7r(s)| = |y - 7r(s)| + \y- 1r(s)\, 

Xi — s — a Xi — s + a 2/! — s — a 2/1 — s + a 



(16) 

|z-7t(s)| ' |x-7ij(s)| |y — ' |2/-7h(s)|* 

necessarily satisfy the following relations: X\ = y\ and x 2 = ±2/2- 

Proof. In order to show this, we will consider (fl"5j) and {jTHJ) as functions of M 3 by 
replacing /i in these expressions with X3 — h. We then transform these expressions 
using the coordinates (|17[) and then set X3 = 2/3 = to prove the lemma. 
Consider the following change of coordinates: 

x\ — s + a cosh p cos 9 y\ = s + a cosh // cos 0' 

(17) X2 = a sinh p sin 8 cos 93 2/2 = ct sinh p' sin 0' cos 93' 

X3 = h + a sinh p sin sin (p y% = h + a sinh p' sin 8' sin 



BISTATIC SAR IMAGING 



9 



where s, a > and h > are fixed and p £ [0, oo), £ [0, 7r] and tp £ [0, 27r). This 
a well-defined coordinate system except for £ = and = 0, 7r. 
In the coordinate system (fl7]l. we have 

| a: — 7t( s )| = a(coshp — cos0), |x — 7i?(s)| = a(coshp+ cos0), 

(J-Oj xi—s — a cosh pcos 6 — 1 x\ — s+a cosh p cos 9+1 

\x— 7t(s) cosh p— cos ' |a:7^(s)| coshp+cosfl 

The terms involving y are obtained similarly. Now (JT5J) and (|16[) transform as 
follows: 

2 cosh p = 2 cosh p 

cosh p cos — 1 cosh p cos + 1 cosh p' cos 0' — 1 cosh p' cos 0' + 1 
cosh p — cos cosh p + cos cosh p' — cos 0' cosh p' + cos 0' 
Using the first equality in the second equation, we have 

cos cos 0' 

cosh 2 p — cos 2 cosh 2 p — cos 2 0' 

This gives cos0 = cos0'. Therefore = 2n7r ± 0', which then gives sin0 = ± sin0'. 
Therefore, in terms of (xi,X2) and (2/1,3/2), we have x\ = y\ and X2 = ±2/2- □ 

Now to finish the proof of the proposition, when x\ = 2/1 and x 2 = 2/2, there is con- 
tribution to WF(K)' contained in the diagonal set A := {(x±, X2, £1, £2; x\, x-x, £1, £2)} 
and when x\ — y± and X2 = —J/2, we have a contribution to WF(K)' contained in 
A, where A := {(x±,X2,£,i,Z,2]x±,-X2,£,i, -£2)}- □ 

We use the following convention in the theorem below. The cotangent variables 
corresponding to x and y are denoted as £ and rj respectively. Then note that 
£ = d x 4> and 77 = —d y 4>. 

We now prove the following theorem: 

Theorem 4.3. Let F be given by ©. Then F*F £ 7 3 >°(A, A). 

Proof. We follow the proof of [BJ Theorem 1.6] closely to prove this result. Recall 
that A and A arc defined by 

A = {xi - 2/1 = x 2 - 2/2 = £1 - f/i = £2 - % = 0}, 
A = {xi - 2/1 = x 2 + 2/2 = £1 - ?7i = £2 + V2 = 0}. 
The ideal of functions that vanish on A U A is generated by 

Px = xi-yi, P2 = x%-yl, P3 = £i-??i, Pa = (x 2 + 2/2) (6 -772), 

= (a?2 -2/2X6 P6=£2- ? ?2- 

Let ^ = gift, for 1 < i < 6, where 51,92 are homogeneous of degree 1 in (£,?/), 
93, 54 and 95 are homogeneous of degree in (£, 77) and % is homogeneous of degree 
— 1 in (£,?/). Let Pi be pseudodiffcrcntial operators with principal symbols pi for 
1 < i < 6. 

In order to prove that F*F £ I p ' l (A,A), for some p, I, we have to show that 
PiK £ H^° c for some so, for 1 < i < 6. By [H Proposition 4.3.1], we have the 
following (up to lower order terms): 

PiK{x,y) = f e i ^ H(s ' v) - R ^ x)) A{x 1 y,s^)p l {x,y 1 ~d x R{s,x)^d y R{ S ,y))A S du 
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We show in the Appendix that each pi can be expressed in the following forms: 

(19) vi = hl ^ V,S) d s <l> + f 12 (x, y, s)d u <t> 

UJ 

(20) p 2 = M^lM^ + /22 ( Xj yi s ) d ^ 

UJ 

(21) p 3 = -d s <t> 

(22) pi = fn(x, y,s)d s (/) + ujf 42 (x,y, s)d u (j> 

(23) p 5 = f 51 (x,y,s)d s (t) + ujf 52 (x,y,s)d u (l) 

(24) p 6 = ujf 61 (x,y,s)d s (f> + uj 2 f 62 (x,y, s)d u <j>. 

Here for 1 < i < 6 and j = 1,2 are smooth functions. 
Therefore 



PxKfay) = J <W**>>rtA(x,v,8,u)) qi ( fu ^ V ' S) d s cf> + f 12 (x, y, s)d u <j> 
= J i).. (> ?±A(x,v,s,u)fu(x,y,8)dsdw 
+ fd u (>^^) ^A(x 7 y,s,uj)f 12 (x,y 7S )d S duj 



dsdcj 



By integration by parts 

jHx,v,s, u ) da (<±A(x, y, s, oj)f n (x, y, s)) dsdu 

\ 1UJ / 



Mi*,v,.*,) du (^A(x, y, s, uj)f 12 (x, y, s) 



dsdcj 



Note that qi is homogeneous of degree 1 in w, and A is a symbol of order 4, hence 
each amplitude term in the sum above is of order 4. 

Therefore by gj Theorem 2.2.1], we have that P X K G Hg c . 

A similar argument works for each of the other five pscudodifferential operators. 
Hence by pH Proposition 1.35], we have that F*F e / P ^(A, A). Because C is a 
local canonical graph away from S, the transverse intersection calculus applies for 
the composition F*F away from S. Hence F*F is of order 3 on A \ S and A\ S. 
Since F*F is of order p + I on A \ £ and is of order p on A \ S, we have that p = 3 
and Z = 0. Therefore the theorem is proved. 

□ 
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Appendix A. 



Here we prove the identities (|T9j) through (|24| that are required in the proof 
of Theorem 14.31 For convenience, and without loss of generality, we will assume 
c = 1. 

In obtaining these identities, it is easier to work in the coordinate system defined 
in (jTT)) . We will work with the extension <fi of the phase function (f> to IR 3 defined 

by 



<t> = — V (2/1 ~ s - a) 2 + y 2 + (j/ 3 - h) 2 + y/( yi -s + a) + y 2 2 + (y 3 - h) 2 

co y v 

- {^J (xi - s - a) 2 + X 2 . + (x 3 - h) 2 + yj {x x - s + a) + x 2 + (x 3 - h) 2 ^ \ . 
Then, using the facts that 

(25) d u <j)\ X3= y 3 =o = d u (j) and d s ^\ X3=V3=0 = d s (f> 7 

we set the third coordinate x 3 = y 3 = to obtain the required identities. 

A.l. Expression for x\ —y\. We now obtain an expression for x\ — y\ of the form 

(26) ^ AiM 8 ^ + /u( ) W 

where /n and /12 are smooth functions. 

That is, denoting A\ = x\ — yx, we would like to obtain an expression of the 
form (|26|) involving d s <f> and 9^0 for 

(27) Ai = a(coshpcos6l — cosh p' cos 0'). 
We have 

dujcj) = 2a(coshp' — coshp) 



and 



L0\ 



( cosh p' cos 9' — 1 cosh p' cos 0' + 1 
\ cosh p' — cos 9' cosh p' + cos 9' 



(28) = 2w 



/ cosh p cos 0—1 cosh p cos + 1 
\ cosh p — cos cosh p + cos 6* 

sinh 2 p cos sinh 2 p' cos 9' 



s cosh p — cos 2 6* cosh p' — cos 2 6*' 
using cosh 2 p cos — cos 9 = sinh 2 p cos 9. Now 



5 

(29) cos(9— ^— = a (coshp' cos 6*' — cosh p cos 9) + a cosh p' (cos 9 — cos 6*') . 
Then 

cos $ — 

(30) A\ = — d^cj) + a cosh/)' (cos 6 — cos 9'). 
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Adding and subtracting c ^ 2 f_ c ™£ 9 , inside ([28]). we have 

sinh 2 p cos 9 sinh 2 p cos 9' 



cosh p — cos 2 9 cosh p — cos 2 9' 
" * ' 



sinh 2 p cos #' sinh 2 p' cos 9' 



cosh 2 p — cos 2 9' cosh 2 p' — cos 2 9' 



ii 

Simplifying / and II. we have, 

2 \f ,i 2 



and 



/ = 



77 



(cos 9 — cos 9') (sinh p)(cosh p + cos 9 cos 6*') 
(cosh 2 p — cos 2 0)(cosh 2 p — cos 2 6>') 

cos 9' sin 2 6*' (cosh p — cosh p')(cosh p + coshp') 
(cosh 2 p' — cos 2 #')(cosh 2 p — cos 2 0') 
cos 9' sin 2 0' (coshp + coshp') d^cj) 
(cosh 2 p' - cos 2 9' ) (cosh 2 p - cos 2 9') 2a 

where we have used the fact sinh 2 p — sinh 2 p' = cosh 2 p — cosh 2 p' . 
Using these calculations, we see 

. . „ „. / (cosh 2 p — cos 2 9) (cosh 2 p — cos 2 9') 

(31) cos 9 - cos 6' = ^ ^- ^ 

\ (cosh p — l)(cosh" p + cos 6* cos 6*') 

(9 S cos 6*' sin 2 0' (coshp + coshp') c^j/m 

~2^T (cosh 2 p' - cos 2 9') (cosh 2 p - cos 2 9') ' ~2a J ' 

Now setting X3 = j/3 = and using fl30[) , we have, 

a cosh p'(cosh 2 p — cos 2 #)(cosh 2 p — cos 2 0') <9 S ^> 

~~ 2" 2" o 

(cosh p— l)(cosh p + cos 9 cos 0') 2a; 

1 / cosh p' cos 0' sin 2 9' (coshp + coshp') (cosh 2 p — cos 2 6*) \ 

— - cos 9 7. 5 • 5 Oujd). 

2\ (cosh 2 p- 1) (cosh 2 p + cos 6> cos 0') (cosh 2 p> - cos 2 9') J 

We can see that no denominator in this expression is zero for 2:3 = (since cosh p > 
1 > cos(? if X3 = 0) and so this expression for A\ is defined and smooth for all 
values of the coordinates. 

We can write A\ in the Cartesian coordinate system as follows. First, for sim- 
plicity, let 



(32) 


Xi = 


\x - 7t(s)| = 


(33) 


x 2 = 


\x-7r(s)\ = 


(34) 


Yi = 


\v-1t{s)\ = 


(35) 


Y 2 = 


\v-1r{s)\ = 



1 (xx - s - 


- a) 2 + x\ ■ 


f h 2 , 


1 (xi - s - 


fa) 2 +i§ ■ 


+ h 2 . 


(yi -s- 


a) 2 + V\ 4 


h 2 , 


(yi - s-f 


■ a) 2 + y\ 4 


-h 2 . 
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Then using these expressions and (fTSf we see that 

rv ( Yi+Y 2 \ ( X 1 X 2 \ ( (X 1 +X 2 ) 2 _ i( Vl -s) 
11 V 2a J \ a 2 I \ ia 2 (Y 1 +Y 2 ) 

'~!+X 2 ) 2 , i(x 1 -s)(y 1 -s) 




\ I3P r- (Xi+x 2 )(yi+y 2 ) 

AlA 2 [^yT) [I - {Yl+ Y 2 )i ) { 



X 1 +X 2 +Y 1 +Y 2 
lot 



Xi+X 2 y :Y .. ( (Xi+*2) 2 A f (X 1+ X 2 ) 2 4 {X1 - S) ( yi - S ) \ 

1112 V 4 " 2 J\ 4q2 ^ (X 1 +X 2 )(Y 1 +Y 2 ) ) 



We see from this second expression for A\ that the coefficient functions are smooth 
in (x,y,s) since the expressions (|32[) - (|35[) are non-zero and smooth. 

A. 2. Expression for x\ — y 2 . Now we write x\ — y\ in the form 

(36) A 2 := x\ - y\ = + /aafo V, 

where /21 and f 22 are smooth functions. A 2 in the coordinate system (|17[) is 

A 2 : = a 2 (sinh 2 p sin 2 6> cos 2 1/3 - sinh 2 p' sin 2 6*' cos 2 ip') 

(37) = a 2 (sinh 2 p sin 2 (9 - sinh 2 p' sin 2 0') 

(38) + a 2 (sinh 2 p' sin 2 0' sin 2 <p' - sinh 2 p sin 2 sin 2 ip). 

For X3 = J/3 = 0, (|38| is 0. So it is enough to obtain an expression of the form (|36| 
for (|37|) . which we still denote by A 2 . 
Using the following identities, 

sinh 2 p = cosh 2 p — 1 and sin 2 = 1 — cos 2 0, 

we have, 

A 2 = a 2 (sinh 2 p sin 2 - sinh 2 p' sin 2 0') 

= a 2 ((cosh 2 p - cosh 2 p') + (cos 2 - cos 2 9') - (cosh 2 p cos 2 - cosh 2 p' cos 2 0')) 
= — — (coshp + coshp')<9 w (/> + a 2 (cos0 + cos 0')(cos — cos0') 

— a (cosh p cos + cosh p' cos 0')^4i 

= — — (coshp + coshp')^^ + a 2 (cos0 + cos 9') (cos 9 — cos0') 

cos 

— a(cosh pcos 6* + cosh p' cos 0')( — d u <fi + a cosh p'(cos — cos 9')) 

= — ^ — (cosh p + cosh p) + (coshp cos 9 + cosh p cos 9') cos 9^j d u (f> 

+ a 2 ^(cos 9 + cos 9') — (cosh p cos 9 + cosh p' cos 0') cosh p'J (cos — cos 0'). 

Now we use the expression for cos0 — cos0' in Equation (|3Tj) and set X3 = 2/3 = 0; 
this shows that x| — y| can be written in the form 

, hi(x,y,s) 

A 2 = d s (j> + f 22 (x, y, s)d w (p. 
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A. 3. Expression for £1 — 771. Now we consider £1 — 771, where we recall that 
(£1,62) = dx4> and (771, T72) = —d y 4> are the cotangent variables corresponding to 
(xi,X2) and (2/1,2/2) respectively. 

Then note that £1 — 771 is + d yi cf>. But this is the same as —d s <fi. Hence 

^3 := £1 - Vi = -d s <t>- 

A. 4. Expression for (#2 — 2/2) (£2 + ^2)- We have (up to a negative sign) 

x 2 x 2 2/2 2/2 



(^2 - 2/2X6 + m) = u{x 2 - 2/2) 
= 2cj 



2a; 



cosh 2 p — cos 2 9 


cosh 2 p' — cos 2 9' 


X2V2 coshp' 


X2V2 coshp 


cosh 2 p' — cos 2 ( 


V cosh 2 p — cos 2 9 


x\ cosh ,0 


x\ cosh p' 



I cosh 2 p — cos 2 6> cosh 2 p 1 — cos 2 0' 

2 _ 2 N^_£°£hp'^ 

cosh p' — cos 2 9' 
X22/2 coshp' £22/2 coshp \ 



x-7t| k-7fll |2/-7t| |y-7ii| 
/ x\ cosh p ?/ 2 cosh p' 

I To 7~9 '. 



cosh p' — cos 2 9' cosh p — cos 2 9 



Here we have added and subtracted cos h ^ p^l'cos 2 e' m * nc P rcv i° us equation. Sim- 
plifying this we get, 

{x 2 - 2/2) (6 + 772) =2uj {x\ - X22/2) 



+ 



(cosh p cosh p' + cos 2 9) (cosh p' — cosh p) 
(cosh 2 p — cos 2 6*) (cosh 2 p' — cos 2 0') 
coshp(cos6' + cos 9') (cos 6* — cos0') 
(cosh 2 p — cos 2 #)(cosh 2 p' — cos 2 9') 



2 2 N COshp' 

+ (x 2 - 2/2)- 



cosh p' — cos 2 9' J 

Now note that cosh p' — cosh p = and we already have expressions for cos 9 - 
cos 9' and x\ — y\ involving combinations of d u (j> and d s (j). 
Hence we can write (X2 — 2/2) (£2 + ^2) in the form 

(X2 -2/2) (6 +%) = hi(x,y, s)d s c/) + ujf42(x,y,s)d cu (j>. 

For future reference, note that our calculation in this section shows that 
(39) 

cosh p cosh p' (cosh p cosh p' + cos 2 6*) (cosh p' — cosh p) 

cosh 2 p — cos 2 6* cosh 2 p' — cos 2 0' (cosh 2 p — cos 2 #)(cosh 2 p' — cos 2 0') 

coshp(cos# + cos 9')(cos9 — cos9') 



(cosh p — cos 2 6>) (cosh pf — cos 2 9') 

A. 5. Expression for (X2 + 2/2X62 — 7 /2)- This is very similar to the derivation of 
the expression we obtained for (x 2 — 2/2) (£2 + ifo)- 
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A. 6. Expression for p\ — rj|- Using ([33]) and (|35|) . we have 

t-2 2 



2 2 _ 2 / / ^2 X 2 \ I V2 2/2 



|»-7t| k-7fl|/ v|y-7rl \v-jr\ 
cosh 2 p 2 cosh 2 p' 



xi 



+ (4 - vl) 



(cosh^ p — cos 2 9) 2 " z (cosh 2 p' — cos 2 6?') 2 

cosh 2 p cosh 2 p' 

(cosh 2 p — cos 2 9) 2 (cosh 2 pi — cos 2 9') 2 

cosh 2 p' 



(cosh p' — cos 2 9') 2 



Now using the computations for x\ — y\ and {x<i — 1/2) (p2 + V2), m particular (|3^|) . 
we can write £| — V2 m ^ ne form 



for smooth functions fei, ft 



62- 



References 

[1] L-E. Andersson, On t/ie determination of a function from spherical averages, SIAM J. Math. 

Anal. 19 (1988), 214-232. 
[2] M. Cheney, A mathematical tutorial on synthetic aperture radar, SIAM Rev., 43 (2001), 301 

- 312. 

[3] M. Cheney and B. Borden, "Fundamentals of radar imaging," Society for Industrial and 

Applied Mathematics (SIAM), Philadelphia, Pennsylvania, 2009. 
[4] J. Duistcrmaat, "Fourier Integral Operators," Birkhauser Boston Inc., Boston, Massachusetts, 

1996. 

[5] J. Cohen and H. Bleistein, Velocity inversion procedure for acoustic waves, Geophysics 44 
(1979), 1077-1085. 

[6] R. Felea, Composition of Fourier integral operators with fold and blowdown singularities, 

Comm. Partial Differential Equations, 30 (2005), 1717-1740. 
[7] R. Felea, Displacement of artefacts in inverse scattering, Inverse Problems, 23, (2007), 1519— 

1531. 

[8] R. Felea and A. Greenleaf, An FIO calculus for marine seismic imaging: folds and crosscaps, 

Comm. P.D.E. 33 (2008), 45-77. 
[9] A. Greenleaf and G. Uhlmann, Non-local inversion formulas for the X-ray transform, Duke 

Math. J. 58 (1989), 205-240. 
[10] A. Greenleaf and G. Uhlmann, Estimates for singular Radon transforms and pseudodifferen- 

tial operators with singular symbols, J. Functional Anal. 89 (1990), 202-232. 
[11] A. Greenleaf and G. Uhlmann, Composition of some singular Fourier integral operators and 

estimates for restricted X-ray transforms, Ann. Inst. Fourier, Grenoble, 40 (1990), 443-466. 
[12] A. Greenleaf and G. Uhlmann, Microlocal techniques in integral geometry, Contemporary 

Math., 113 (1990), 121-136. 
[13] V. Guillcmin, Cosmology in (2 + 1)- dimensions, cyclic models, and deformations of Mi,\, 

Annals of Mathematics Studies, 121. Princeton University Press, Princeton, NJ, 1989. 
[14] V. Guillcmin and S. Sternberg, Geometric Asymptotics, American Mathematical Society, 

Providence, RI, 1977. 

[15] V. Guillcmin and G. Uhlmann, Oscillatory integrals with singular symbols Duke Math. J. 
1981, 48 (1), 251-267. 

[16] L. Hormandcr, Fourier integral operators, I, Acta Mathematica 127 (1971), 79-183. 
[17] A. M. Home and G. Yates, Bistatic synthetic aperture radar Proc. IEEE Radar Conf., 2002, 
6-10. 

[18] L. Hormander, The Analysis of Linear Partial Differential Operators, I -IV, Springer Verlag, 
New York, 1983. 



16 



VENKY P. KRISHNAN AND ERIC TODD QUINTO 



[19] A.K. Louis and E.T. Quinto, Local Tomographic Methods in SONAR, Surveys on Solution 

Methods for Inverse Problems (Vienna/New York) (D.Colton, H.Engl. A.Louis, J. McLaughlin, 

and W. Rundell, eds.), Springer Verlag, 2000, pp. 147-154. 
[20] R. Melrose, G. Uhlmann, Lagrangian intersection and the Cauchy problem Comm. Pure Appl. 

Math., 1979, 32 (4), 483-519. 
[21] C. J. Nolan and M. Cheney, Synthetic aperture inversion, Inverse Problems, 18, 2002, 221- 

235. 

[22] C. J. Nolan and M. Cheney, Microlocal analysis of synthetic aperture radar imaging, J. Fourier 

Anal. Appl., 10, 2004, 133-148. 
[23] E.T. Quinto, The dependence of the generalized Radon transform on defining measures, Trans. 

Amcr. Math. Soc. 257 (1980), 331-346. 
[24] F. Treves, "Introduction to pseudodiffercntial and Fourier integral operators. Vol. 1 and 2," 

Plenum Press New York 1980. 
[25] C. E. Yarman, B. Yazici and M. Cheney, Bistatic synthetic aperture radar imaging with 

arbitrary trajectories, IEEE Transactions on Image Processing, 17, (2008), 84-93. 

E-mail address: venkyp.krishnan@gmail.com 
E-mail address: todd.quinto@tufts.edu 



